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The adaptive immune system reacts against pathogenic nonself, whereas it normally 
remains tolerant to self. The initiation of an immune response requires a critical antigen(Ag)- 
stimulation and a critical number of Ag-specific T cells. Autoreactive T cells are not 
completely deleted by thymic selection and partially present in the periphery of healthy 
individuals that respond in certain physiological conditions. A number of experimental and 
theoretical models are based on the concept that structural differences discriminate self 
from nonself. In this article, we establish a mathematical model for immune activation in 
which self and nonself are not distinguished. The model considers the dynamic interplay 
of conventional T cells, regulatory T cells (Tregs), and IL2 molecules and shows that the 
renewal rate ratio of resting Tregs to naive T cells as well as the proliferation rate of acti- 
vated T cells determine the probability of immune stimulation. The actual initiation of an 
immune response, however, relies on the absolute renewal rate of naive T cells. This result 
suggests that thymic selection reduces the probability of autoimmunity by increasing the 
Ag-stimulation threshold of self reaction which is established by selection of a low num- 
ber of low-avidity autoreactiveT cells balanced with a proper number of Tregs. The stability 
analysis of the ordinary differential equation model reveals three different possible immune 
reactions depending on critical levels of Ag-stimulation: a subcritical stimulation, a thresh- 
old stimulation inducing a proper immune response, and an overcritical stimulation leading 
to chronic co-existence of Ag and immune activity. The model exhibits oscillatory solutions 
in the case of persistent but moderate Ag-stimulation, while the system returns to the 
homeostatic state upon Ag clearance. In this unifying concept, self and nonself appear 
as a result of shifted Ag-stimulation thresholds which delineate these three regimes of 
immune activation. 
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INTRODUCTION 

The immune system is continuously exposed to a wide variety of 
disturbances. Such disturbances are recognized by T cells via anti- 
gen presentation. Antigen presentation is a process in which anti- 
gen presenting cells (APC) capture the antigens, break them into 
small peptides, couple them with MHC molecules, and present 
them on the cell surface, thus enabling their recognition by T cells 
( 1-3 ) . The majority of disturbances that the immune system deals 
with are pathogenic nonself-antigens. Since the APCs break down 
the nonself-antigens into smaller peptides and present them on 
their surface, presented peptide of nonself might have overlaps 
with self-peptides (4, 5). 

In addition, rapidly evolving nonself pathogens, such as Hepati- 
tis C virus, might acquire similarities to self-antigens (6). Apart 
from nonself, altered self such as cancer cells is also a distur- 
bance that has to be recognized by the immune system. Therefore, 



an ideal immune system has to find a solution for dealing with 
nonself, self-similar nonself, and self-disturbances (Figure 1). 

As a general solution, the immune system generates T cell clones 
with random specificities that could potentially recognize any pep- 
tides, including self-peptides. The classical idea that the T cell 
repertoire has to be self-tolerant and T cells should not react to 
self-peptides, assumes that self-reactive T cells should be elimi- 
nated. This assumption is partially true, as T cell clones which fully 
recognize self-peptides in the thymus undergo clonal deletion, in 
the so-called negative selection process (7, 8). 

The self-tolerance resulting from negative selection is called 
central tolerance. A stringent central tolerance induction and dele- 
tion of all autoreactive T cells is believed to create holes in the 
specificity space of the T cell repertoire (9, 10) by prohibiting 
immune responses against self-similar nonself and altered self. 
Hence, a too stringent central tolerance does not seem beneficial. 
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FIGURE 1 | Conceptual figure of different disturbances in the immune 

system. Pathogenic nonself-disturbances are recognized and attenuated by 
nonself-specific T cells. However, recognition and attenuation of altered self 
and self-similar nonself-disturbances is challenging for the immune system 
due to the existence of self-tolerance mechanisms; without self-specific 
immune cells, the immune system is not able to initiate an immune 
response against these disturbances. 



In line with this idea, there is evidence that negative selection only 
partially deletes autoreactive T cells because availability of self- 
peptides required for negative selection in the thymus is limited 
and T cells spend only a limited time in the process of thymic 
selection (11-13). Autoreactive T cells escaping negative selection 
have been shown to be involved in autoimmunity (14). They nor- 
mally exist in healthy individuals and are quiescent in steady state 
conditions in the presence of their cognate self- antigen (15). 

Escaped autoreactive T cells are under the control of peripheral 
tolerance. A prominent mechanism of peripheral tolerance among 
others [reviewed in Ref. (16)] is induced by CD4+ Foxp3+ regu- 
latory T cells (Tregs) (17). The majority of these cells, known as 
natural Tregs, are hypothesized to be selected from autoreactive T 
cells in thymus (18, 19). The main role of Tregs is the regulation 
of the immune response by suppression of the effector functions 
of conventional T cells (Tconv). 

Despite the necessity of suppression by Tregs for avoiding 
autoimmunity (20, 21), production of too large numbers of Tregs 
in the thymus might prevent beneficial effector responses. There- 
fore, a too stringent peripheral tolerance induction by selection of 
large numbers of Tregs in the thymus does not seem favorable. 

In view of this background, how does the immune system 
balance the tolerance mechanisms in order to ensure immune 
responses to any kind of disturbances including self-disturbances, 
yet staying tolerant to self in healthy homeostasis? Here, we address 
this question by using a mathematical model of immune activa- 
tion that relies on identical components for self and nonself, i.e., 
using the same set of ordinary differential equations. The model 
considers the thymic production of Tregs and Tconvs as well as 
the dynamic interplay between Tregs, Tconvs, and IL-2 molecules 
in the presence of antigen(Ag)-stimulation in the periphery. The 
model is exploited to reveal the parametric regime of the immune 
system in which an immune response against self is restricted, but 
not impossible. 

The interplay between Tregs and Tconvs during immune 
responses is a topic of extensive mathematical modeling (22-28). 
Leon and co-workers (22) proposed a series of models for studying 
immune tolerance by considering APCs, Tconvs, and Tregs. Their 
models rely on the assumption that regulatory interaction between 



Tregs and Tconvs takes place only in simultaneous conjugation 
with an APC. As a result of this assumption, efficient suppression 
of Tconvs requires a minimum population of Tregs per APC (29). 
As an extension, a crossregulation model is proposed by Carneiro 
and co-workers (26) in an attempt to incorporate Tregs in a coher- 
ent theory of the immune system. According to their model that 
shows a bistable behavior, immunity to a given Ag arises as compet- 
itive exclusion of Tregs by the expansion of Tconvs and tolerance 
results from limited APC availability or above threshold Treg num- 
bers. Since the interactions between Tregs and Tconvs is assumed 
to depend on the density of the APCs, increasing the APC avail- 
ability decreases the simultaneous conjugate formation of Tregs 
and Tconvs with the same APCs and hence, it is sufficient to break 
the immune tolerance. 

An alternative concept was brought forward in a model pro- 
posed by Carneiro and co-workers (23) that assumes APC- 
independent interactions between Tconvs and Tregs for immune 
suppression which will be also used in our model. A direct inter- 
action of Tconvs and Tregs was identified by experiments (30). 
The authors concluded that efficient immune suppression still 
requires a minimum population of Tregs regardless of the number 
of APCs. 

In contrast to the aforementioned studies, we do not consider 
the conjugate formation of Tregs and Tconvs with APCs and there- 
fore, there is not a competition between these cells for Ag. Instead, 
the role of APCs is indirectly captured by an Ag-stimulation factor 
which is the activation rate of naive T cells and resting Tregs with 
identical Ag-specificity by APCs bearing their cognate Ag. In addi- 
tion, we explicitly consider the dependency of Tregs on Tconvs 
through the growth factor IL-2. 

Burroughs and co-workers (24) investigated Treg-induced inhi- 
bition of cytokine secretion by effector T cells. By assuming that 
Tregs are activated by self Ag and locally maintained by nonlinear 
competition for tissue-derived cytokines that are solely utilized 
by Tregs, the authors analyzed the role of local active Treg popula- 
tion size in the balance between suppressor and effector responses. 
Stimulation of Tregs and Tconvs is described by independent para- 
meters. In contrast to their model, thymic output maintains the 
homeostatic population of Tregs in our model. Another essential 
difference is that Ag-stimulation of Tregs and Tconvs is described 
with a unified self-nonself concept and Tregs are assumed to also 
respond to nonself Ag-stimulation (31). 

Parametric steady state analysis of the model provides insights 
about the physiological range of model parameters, and deter- 
mines the overall conditions under which immune responses 
against self are possible. Furthermore, the impact of model para- 
meters on the requirements for the initiation of immune reactions 
against self is analyzed. The model proposes that disturbed home- 
ostatic balance between autoreactive T cells and Tregs increases 
the susceptibility to autoimmunity or cancer. 

RESULTS 

The mathematical model is constructed starting from a simple 
model of the immune response including essential components 
only. Then, additional complexity is incrementally added to the 
model to a degree allowing for validation and analysis of tol- 
erance versus immunity. The scheme of the complete model is 
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FIGURE 2 | Model of dynamic interplay between conventional T cells and 
regulatory T cells. Nonself-specific as well as some self-specific thymocytes 
that survived negative selection and were not selected asTregs enter the 
periphery as na'ive T cells. A part of detected autoreactive thymocytes 
differentiate intoTregs in the thymus and reside in the periphery in resting 
state. Upon Ag-stimulation by APCs, na'ive T cells and resting Tregs become 



activated. In contrast to activated T cells, activated Tregs do not secrete IL-2, 
but both activated populations proliferate in dependence on the presence of 
IL-2 (46). Activated Tregs suppress activated T cells in a cell-contact- 
dependent and cytokine-driven manner. Activated T cells undergo Fas- 
induced apoptosis by interacting with each other (fratricide). In contrast, Tregs 
are resistant to Fas-induced apoptosis (68). 



depicted in Figure 2. The model is conceptually independent of 
the self/nonself nature of the immune response, and differences of 
the immune responses against self versus nonself are reflected in 
different parameter values of the same model. 

AN IMMUNE RESPONSE REQUIRES SUFFICIENT DIVISION AND IL-2 
SECRETION RATE OF ACTIVATED T CELLS 

Immune responses arise from massive proliferation of activated T 
cells and their subsequent effector function. Our simplest model 
attempts to capture the dynamic characteristics of an activated T 
cell population (T): 

fdT 

— = alT- bT 

dt (1) 
dl 

— = dT - elT - fl 
I dt 

Activated T cells have a mean lifespan lib and secrete IL-2 (7) 
with rate d. Available IL-2 decays with rate / and is consumed by 
activated T cells with rate e. Activated T cells are able to proliferate 
(with rate al) in the presence of IL-2. This IL-2 dependent prolif- 
eration rate is considered as a linear function of IL-2 in model (1). 
The impact of considering a nonlinear proliferation rate (a Hill- 
function of IL-2) instead of the linear term alT is given in Section 
"Nonlinear Proliferation Rate of Conventional and Regulatory T 
Cells" in Appendix. 

Steady state analysis of the model ( 1 ) is given in Section "Steady 
State Analysis of Model (1)" in Appendix. This model has two 
equilibrium points: 

(Ii.Ii) = (0,0), (T 2 ,I 2 ) =( hf , (2) 
\ ad — be a J 

By assuming the biological range of parameters (all parameters 
are positive), the trivial equilibrium point ( T\ , I\ ) is stable and the 



nontrivial equilibrium (T2, 12) is unstable. T 2 is positive if and 
only if: 

ad-be > 0 (3) 

The unstable equilibrium point imposes a threshold for initial 
conditions of the model in which the activated T cells prolifer- 
ate unlimitedly, which in this simplest model, corresponds to an 
efficient immune response. This can be visualized by the phase 
portrait of the model as shown in Figure 3A. The condition (3) 
imposes a quality constraint on activated T cell clones to enter a 
highly proliferative state and implies that among T cell clones that 
are in the activated state, only the T cell clones with a sufficiently 
high proliferation rate (a) or IL-2 secretion rate (d) are able to 
contribute to the immune response against Ag. Since both, the 
proliferation and IL-2 secretion rate of activated T cells depend on 
the affinity/ avidity of their TCR to the presented Ag (32-34), con- 
dition ( 3 ) implies that the existence of T cell clones with sufficiently 
high specificity for the presented Ag is required for induction of 
an immune response. Similar implications were derived from a 
model that considers a nonlinear IL-2 dependent proliferation rate 
of activated T cells (Nonlinear Proliferation Rate of Conventional 
and Regulatory T Cells in Appendix). 

The major focus of central tolerance is to eliminate T cells that 
are self-specific. Therefore, it is unlikely that highly self-specific 
T cells escape from central tolerance, as they are more effectively 
detected and eliminated in the thymus (12, 34). It is thus expected 
that autoreactive T cells in the periphery are less aggressive than 
the ones that undergo clonal deletion in the thymus, and may not 
fulfill condition (3). 

INITIATION OF AN IMMUNE RESPONSE REQUIRES A MINIMUM 
HOMEOSTATIC POPULATION OF NAIVE T CELLS AND ANTIGEN 
STIMULATION 

Continuous thymic production of na'ive T cells maintains the 
peripheral number and diversity of mature naive T cells (35), 
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although other mechanisms such as stimulation of T cells with 
self-antigens and IL-7 have been shown to be involved (36). Upon 
Ag-stimulation by activated APCs, naive T cells with high avidity to 
the presented Ag become activated. Here, we take into account the 
dynamics of the naive T cell population (N) and T cell activation 
by Ag-stimulation (P), as described in equations (4). We assume 
that naive T cells with identical Ag-specificity have a homeosta- 
tic population in the periphery that is established by naive T cell 
renewal (by rate No) and natural cell death (with rate g): 



dN 
~dT 
dT 
It 1 
dl 
dt = 



= /(N) = N 0 -gN-/3N 
= alT -bT + pN 

dT - elT - fl 

T cell activation k(t) is defined as 

k(t) = pN{t) 



(4) 



(5) 



Steady state analysis of model (4) is given in Section "Steady 
State Analysis of Model (4)" in Appendix. This model has either 
2 or no equilibrium points dependent on the steady state value of 
T cell activation (k). According to the bifurcation diagram of the 
model depicted in Figure 3B, which is obtained by treating k as 
bifurcation parameter, model (4) has no equilibrium points for: 



k > le- 



ad/ 



be 
ad 



(6) 



which corresponds to the unlimited proliferation state of activated 
T cells. Therefore, condition (6) has to be satisfied for initiation of 
an immune response. However, according to model (4), the steady 
state value of T cell activation ( k) is limited by naive T cell renewal 
(No) and Ag-stimulation (ft): 



g + p 



(7) 



Therefore, according to equations (6) and (7), there exists an 
Ag-stimulation range 



gk- 



N 0 - k_ 

in which an immune response is initiated if: 

No > k_ 



(8) 



(9) 



Condition (9) implies that the renewal rate of naive T cells 
plays a critical role for the initiation of immune responses. In 
other words, without a sufficient renewal rate of naive T cells, 
the immune response cannot be initiated by any Ag-stimulation. 
Instead, Ag-stimulation results in a subcritical immune response 
which is interpreted as insufficient for pathogen clearance. By 
increasing the proliferation rate or IL-2 secretion of activated T 
cells or the renewal rate of naive T cells, the threshold of Ag- 
stimulation required for initiation of an immune response is 
decreased [equations (6) and (8)]. Therefore, central tolerance 
is able to tune the initiation criterion of self reaction not only 
by limiting the quality of autoreactive T cells, but additionally 
by restricting the renewal rate of autoreactive T cells. As central 
tolerance does not limit nonself-specific T cells, according to the 
model, these cells exhibit a lower threshold of activation by nonself 
Ag-stimulation. 

FRATRICIDE: A MECHANISM TO LIMIT BUT NOT TO SUPPRESS 
IMMUNE RESPONSES 

The immune response in model (4) is characterized by unlimited 
proliferation of activated T cells which is physiologically unreal- 
istic. The linear death term of natural death of activated T cells 
in model (4) is not sufficient to limit proliferation, and requires a 
nonlinear limiting factor. A potential mechanism of limiting the 
immune response is activation-induced cell death (AICD) in acti- 
vated T cells, known as fratricide (37). Upon T cell activation, 
death ligand (FasL) and receptor (Fas) proteins are expressed on 
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FIGURE 3 | (A) Qualitative phase portrait of model (1): the stable manifold of 
saddle node defines a threshold for the initial conditions that allow for 
unlimited proliferation of activated T cells. (B) Bifurcation diagram of 



model (4) by treating k as bifurcation parameter. Stable and unstable 
equilibrium points are shown by black and red lines, respectively. For k > k_, 
the immune response enters the regime of unlimited proliferation. 
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the surface of T cells. Followed by expression of these proteins, 
T cells eliminate themselves in a cell-contact-dependent manner. 
The fratricide mechanism is modeled by a nonlinear death term 
(cT 2 ), as proposed by Callard et al. (37): 



dN 

~dt 

dT 

It 

dl 

dt = 



= f(N) =N 0 -gN-/3N 
-- alT -bT- cT 2 + pN 
dT - elT - fl 



(10) 



The steady state analysis of model (10) is provided in Section 
"Steady State Analysis of Model (10)" in Appendix. This model 
has either 3 or 1 equilibrium points, depending on the value of 
fratricide death rate c. The bifurcation diagram of the model (10) 
with respect to c is depicted in Figure 4A for (p = 0). When c 
satisfies 



(11) 



f- 1 (Vad - */be\ 



the stable equilibrium point (73) exists and corresponds to a sat- 
urated population of activated T cell. When the conditions (3) 
and (11) are fulfilled, the model (10) exhibits the bifurcation dia- 
gram plotted in Figure 4B with respect to the steady state value 
of T cell activation (k). The fratricide mechanism added a large 
stable equilibrium point (T3) to the model which imposes a sat- 
uration level to the activated T cell population. The larger the c, 
the smaller the saturated population of activated T cells is. Simi- 
lar to model (4), model (10) shows an initiation threshold of the 
immune response (k > kf). Despite solving the issue of unlim- 
ited proliferation of activated T cells by the fratricide mechanism, 
model (10) bears a hysteresis characteristic so that the immune 



H 
H 

4 

a 10 



10 10 

Fratricide death rate (c) 



response cannot be suppressed when Ag-stimulation (P) is 
removed. 



DYNAMIC INTERPLAY OF ACTIVATED T CELLS AND TREGS 

Tregs are essential in maintaining self-tolerance and immune 
homeostasis by preventing autoimmunity and limiting chronic 
inflammation in the periphery. However, they might also limit 
beneficial responses by inducing tolerance to pathogens (38, 39) 
or limiting anti-tumor immunity (40, 41). One functional role of 
Tregs is to shut down the cell-mediated immune response via cell- 
contact-dependent and inhibitory cytokine-driven suppression of 
activated T cells (42). Two different subsets of Tregs were iden- 
tified. Natural Tregs are the dominant subset of peripheral Tregs 
(43) and are selected in the thymus. In our model, we consider 
only natural Tregs and neglect the induced Treg subset that dif- 
ferentiates from naive T cells. Like for naive T cells, the thymus 
contributes to the renewal of resting Tregs (N) by continuously 
selecting them from thymocytes. The renewal of resting Tregs is 
assumed to occur by rate No. Since we are interested in the relative 
renewal of resting Tregs and naive T cells, we assume that: 



N 0 = A No 



(12) 



Tregs remain in the resting state until they are stimulated by 
Ag (p) and become activated in a TCR-dependent manner. The 
dynamic population of the resting Treg compartment is assumed 
to be the same as the naive T cell compartment in (4) and (10) 
(dN/dt = /(N)). Activated Tregs (_R) are assumed to suppress 
activated T cells (by rate y ) . Survival and proliferation of activated 
Tregs depends strictly on IL-2, produced by activated non-Tregs 
(44-46). The IL-2 dependent proliferation rate of Tregs is con- 
sidered as a linear function of IL-2 (see Nonlinear Proliferation 
Rate of Conventional and Regulatory T Cells in Appendix for a 
nonlinear case). In contrast to activated T cells, activated Tregs 
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FIGURE 4 | (A) Bifurcation diagram of model (10) with p = 0 using the 
fratricide death rate c as bifurcation parameter. No immune response 
exists for fratricide death rates larger than c_ due to extensive 
activation-induced cell death. The trivial equilibrium point is omitted in 
this figure. (B) Bifurcation diagram of model (10) using k as bifurcation 



parameter: an immune response can be initiated for large values of k. 
However, due to hysteresis characteristic in this model, the immune 
response is not suppressed after decreasing T cell activation (k). Stable 
and unstable equilibrium points are shown by black and red lines, 
respectively. 
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lack the ability to secrete IL-2 (47). The relative proliferation 
rate of activated Tregs and activated T cells is controlled by the 
parameter e: 



dT 
dt 
dR 
dt 
dl 
dt 



= alT -bT - cT 2 - yRT + @N 

= eaIR- bR + 0 (13) 

= dT - eI(T + R) - fl 



The parameters are given in Table 1 and the model components 
are illustrated in Figure 2. Treg activation k(t) is defined as 



k(t) = pN(t) 



(14) 



According to equations (7) and (12), the steady state value of 
Treg activation (k) is given by 



PN 0 



kk 



(15) 



The equilibrium points of model (13) are given in Section 
"Steady State Analysis of Model (13)" in Appendix. By incor- 
porating the Treg compartment to model (10), two additional 
equilibrium points (T4 and T5) emerged for p = 0. The equilib- 
rium point of interest (T4), which depends on the Treg-associated 
parameters (e, y), has an impact on the topological changes of 
the phase portraits of the model under variations of the bifur- 
cation parameter k. The value of e and y are assumed to be in 
a range where the model does not inherit the hysteresis char- 
acteristics of immune responses from model (10) in which the 
immune response is not suppressed after resolving Ag-stimulation 
(P ) . Then, the bifurcation diagrams of model ( 1 3 ) for two different 
values of k are obtained by treating k as the bifurcation parameter 
(Figure 5). Depending on the value of k, the model has either 5 or 
3 equilibrium points. 

By varying the relative renewal rates of resting Tregs and naive 
T cells [k in equation (12)] a k t ^ can be found, so that no immune 



response can be initiated for any value of k, if k > k t ^ (Figure 5A). 
For k < k t h (Figure 5B), there exists a T cell activation threshold 
(fc,) such that for k > k\ the immune response can be initiated. 
However, in contrast to model (10), the immune response is com- 
pletely suppressed by activated Tregs if k decreases to a lower value 
than kj (gray region in Figure 5B). For persistent Ag-stimulation 
with k > kj, two scenarios are possible. An oscillating immune 
response is induced when k remains in the range of k, < k < k s 
(red region in Figure 5B). For k> k s the immune response is sup- 
pressed after its initiation to a minor immune response with an 
activated T cell population T4 due to over-suppression of activated 
T cells by over-activation of Tregs (magenta region in Figure 5B). 
In the latter case (k > k s ), despite proper T cell stimulation, only a 
minor immune response is induced (and antigen is not cleared). 
Instead a chronic co-existence of antigen and inefficient immune 
activity is established. Therefore, according to the model, a range of 
T cell and Treg activation (ki<k < k s ) exists in which an efficient 
immune response is induced. Outside of this range, the antigen 
persists because of under-stimulation of naive T cells, or over- 
stimulation of Tregs. According to equation (7), the existence of 
Ag-stimulation thresholds p, and p s which correspond to the val- 
ues of ki and k s , respectively, depends on the renewal rate of naive 
T cells (No); Pi exists if No > h and p s exists if No > k s . Increas- 
ing the renewal rate of naive T cells reduces the Ag-stimulation 
required for initiation( / S 1 )/over-suppression(/J s ) of the immune 
response. 

The peak immune response depends on the value of the 
Treg-associated equilibrium point (T4) which in turn depends 
on Treg-associated parameters. However, the fratricide-associated 
equilibrium point (T$) is a limiting factor for the maximum 
population of activated T cells if the fratricide death rate (c) is 
sufficiently high. 

According to our model, sufficient activated Tregs are required 
to suppress the proliferative response of activated T cells. These 
are supplied by two processes: Treg activation (k) which depends 
on Ag-stimulation (/?), and Treg proliferation which depends on 
the IL-2 secretion by activated T cells. With a low Ag-stimulation 
and insufficient Treg activation (k = pN), Treg proliferation has 
to account for immune suppression. Since Treg proliferation is 



Table 1 | Parameters used for model analysis. 



Parameter Value Description Dimension 



a 


0.4 


Proliferation rate of activated T cells 


molecules -1 time -1 


b 


0.1 


Natural death rate of activated T cells and Tregs 


mime -1 


c 


10" 5 


Fratricide death rate of activated T cells 


cells -1 time -1 


d 


0.01 


IL-2 secretion rate by activated T cells 


molecules cells -1 time 


e 


0.01 


IL-2 consumption rate by activatedT cells andTregs 


cells -1 time -1 


f 


1 


IL2 decay rate 


time -1 


9 


B 


Natural death rate of naive T cells and resting Tregs 


time -1 


P 


[0,oo) 


Ag-stimulation of naive T cells and resting Tregs 


time -1 


y 


0.1 


Treg-mediated suppression rate 


cells -1 time -1 


€ 


0.6 


Proliferation rate ratio Treg/Tconv 




N 0 


4 


Renewal rate of naive T cells 


cells time -1 


I 


0.006, 0.02 


Relative renewal rate of resting Tregs and nai'veT cells No/ No 




No 


X No 


Renewal rate of resting Tregs 


cells time -1 
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FIGURE 5 | Bifurcation diagram of model (13) using k as the 
bifurcation parameter with (A) X = 0.02 and (B) X = 0.006 Stable and 
unstable equilibrium points are drawn by black and red solid lines, 
respectively. Dashed black lines represent the stable limit cycles by 
showing the maximum and minimum populations of oscillating activated T 
cells for persistent k. Depending on the values of X and k, an immune 



dependent on the availability of IL-2, sufficient activated T cells are 
required to secrete IL-2 and induce immune suppression. There- 
fore, activated T cells undergo the proliferation up to a level that 
sufficient IL-2 is available for Treg proliferation and subsequent 
immune suppression. In contrast, by facilitated Treg activation 
(fc), less Treg proliferation is required for suppressing activated 
T cells which means that the dependency of immune suppres- 
sion on proliferation of activated T cells decreases. Consequently, 
by increasing Ag-stimulation (/S) in the range of /3, < ft < /3 S 
(red region in Figure 5B), Treg activation (k) increases as well 
which results in a reduced maximum population of activated T 



response is not initiated (gray), is initiated (red) or over-suppressed 
(magenta). With parameter values given in Table 1, the threshold becomes 
= 0.01183. The time-courses of the activated T cell population 7~(t) were 
deduced from a numerical solution of model (13) with zero initial 
conditions and persistent p. The unstable negative equilibrium point (7 5 ) is 
not shown in the plots. 



cells (Figure 5B, dashed black line) and an increased frequency 
of oscillations. By further increasing Ag-stimulation to ft > fi s 
(magenta region in Figure 5B), Treg activation (fc) completely 
prevent oscillating immune response. 

In the same way, by increasing the relative homeostatic popu- 
lation of resting Tregs and naive T cells (A. > A.^), Treg activation 
increases up to a level that Treg suppression does not depend on the 
proliferative response of activated T cells. Thus, activated T cells are 
not able to enter the massive proliferation for any Ag-stimulation 
level, as shown in Figure 5A. Similar results were derived from a 
model that considers a nonlinear IL-2 dependent proliferation rate 
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of activated T cells and Tregs (see Nonlinear Proliferation Rate of 
Conventional and Regulatory T Cells in Appendix). 

DISCUSSION 

In this study, a model of the dynamic interplay between effec- 
tor and regulatory immune responses was examined to investigate 
the requirements for the initiation of an immune response by 
Ag-stimulation. The model unifies several components developed 
in previous studies, such as IL-2 dependent proliferation of T 
cells (48), fratricide-induce programed cell death (37), IL-2 com- 
petition between activated T cell and activated Tregs (24), and 
Treg-mediated immune suppression (23, 24, 28). Homeostatic 
division of T cell compartments was not considered in the present 
study, such that the main renewal source of T cells in the absence 
of Ag-stimulation is the thymus. While the presented model is 
still simplifying the real situation in many aspects, the stability 
analysis revealed a number of reasonable results matching many 
experimental findings and allowing for an analysis of reasons for 
immune failure and autoimmunity. 

The model predicts three qualitatively different immune 
responses depending on the level of antigenic stimulation. At first, 
a threshold stimulation p, is required in order to get an immune 
response at all. Secondly, in a limited range of Ag-stimulation 
P € (Pi, p s ) an efficient immune response is induced. Tregs limit 
the duration of the immune response. If the antigen was cleared by 
the first immune response, further immune activity would be sup- 
pressed by Tregs. However, if the first peak of the immune response 
fails to clear the antigen, but keeps the antigen in the stimulation 
range /3,- < p < p s , the immune system attempts to clear the anti- 
gen with subsequent immune responses, which corresponds to the 
oscillatory solution depicted in Figure 5B. If the immune response 
failed to control the antigen spread, antigenic stimulation would 
be further increased to p > p s , leading to the third class of immune 
responses. Tregs are over-stimulated and suppress immune activ- 
ity. In this situation, a chronic persistence of the antigen would 
develop. Treg-mediated over-suppression of immune responses 
in chronic infections is well-established (reviewed in Ref. (49)). 
According to our model, depletion of resting Tregs restores the 
immune response by transiently decreasing X and by this increas- 
ing p s . This notion is consistent with the experimental model of 
chronic infections according to which depletion of Tregs results 
in the restoration of effector immune response and restriction 
of antigen spread (50, 51). A key feature of our model is that 
the immune response does not rely on a stable equilibrium point 
with a dominant population of activated T cells which is typi- 
cally derived from existing bistable models. It rather relies on a 
transient response (or stable limit cycles in the case of persistent 
Ag-stimulation) which originates from T-cell-mediated suppres- 
sion and IL-2 consumption by Tregs. Moreover, the role of Tregs 
in the chronic state of the immune response is not represented by 
available models. 

According to our model, the qualitatively different immune 
responses and their requirements are dependent on the quality 
and quantity of Tconv and Treg clones responding to the Ag- 
stimulation. The proliferation rate of activated T cells, which 
depends on their avidity to the stimulating antigen determines the 
existence of an Ag-stimulation threshold (Pi) which is required 



for the initiation of an immune response. The absolute renewal 
rate of naive T cells (No) adjusts the Ag-stimulation threshold 
which exists when the renewal rate of resting Tregs remains below 
a threshold value (X < X t h)- Further Treg-associated parameters, 
namely the proliferation rate of Tregs (e ) and the Treg-mediated 
suppression rate (y), also affect the existence and the level of the 
Ag-stimulation required for initiation (pi) and over-suppression 
(Ps) of immune responses. By increasing the proliferative (e) and 
suppressive (y) activity of Tregs, p, increases, whereas p s decreases 
up to a level where the initiation of an immune response is 
completely impossible for any Ag-stimulation. Interestingly, when 
proliferation rate of activated Tregs exceeds the one for activated 
T cells (e > 1) a massive proliferation of activated T cells is still 
required for subsequent immune suppression by Tregs. Thus, IL- 
2 secretion by sufficiently large numbers of activated T cells is a 
strict requirement for immune suppression. Note also that with- 
out Tregs, a return to the homeostatic state is not possible, even 
when the antigen was cleared. 

Considering all aforementioned parameters controlling the ini- 
tiation of an immune response, is it beneficial for the immune sys- 
tem to completely avoid self reaction, or is there a benefit in allow- 
ing self reaction? Clearly, autoreactive T cells exist in the periphery 
of healthy individuals as a normal component of the T cell reper- 
toire (12, 14, 52, 53). These cells respond to self-tissue destruction 
even in the presence of Tregs and without genetic predisposition to 
autoimmunity (15). Although the activation of autoreactive T cells 
has been shown to be involved in autoimmunity (12), several lines 
of evidences indicate that these cells are required for limiting self- 
destruction by supporting self-regenerative processes (54-56). In 
addition, the anti-tumor immune responses evoked by autoreac- 
tive T cells are beneficial (34, 57). Therefore, it seems unlikely that 
autoreactive T cells escaping from the thymus are simply a result of 
thymic selection error that can disturb self- tolerance under certain 
physiological conditions. Instead, these evidences imply that ben- 
eficial self reaction is allowed in the immune system. According to 
the mathematical model, immune reactions against self are only 
possible with a critical homeostatic population of autoreactive T 
cells (or sufficient renewal rate No) which is balanced by a proper 
number of Tregs (X < A t h) which corresponds to region (C) or 
(D) in Figure 6. Since the T cell repertoire is normally stimulated 
with an endogenous level of self-antigens in the periphery which 
does not evoke any self reaction, the Ag-stimulation threshold for 
initiating an immune response (pi) should be sufficiently high in 
comparison to a typical nonself Ag-stimulation. According to our 
model, this is achieved by ensuring a low renewal rate (No ) of low- 
avidity autoreactive T cells and a high, but balanced renewal rate 
of Tregs (high X but lower than X t h). In other words, according to 
Figure 6, by choosing No close to k, and higher value of which 
is obtained by higher X, a large Ag-stimulation threshold (Pi) for 
the initiation of immunity against self can be achieved. 

Aging of the immune system, the so-called immunosenescence, 
is characterized by involution of thymus, decreased number of 
thymic output, contraction in T cell diversity, and disturbed T cell 
homeostasis which all result in attenuated immune function and 
susceptibility to infectious diseases and cancer in the elderly (58, 
59). By decreasing thymic output, the homeostatic population of 
some T cell clones diminishes which leads to the creation of holes 
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renewal rate of naive T cells N 0 



FIGURE 6 | The balance between renewal rate of naive T cells and 
resting Tregs. The relative renewal rate of na'ive Tcells and restingTregs 
(X = Wo/Wo) determines the existence of an immune response. The 
initiation of an immune response requires a sufficient renewal rate of 
naive T cells (W 0 ). (A) For X > X th , the immune response does not exist for 
any value of W 0 and Ag-stimulation (j8). (B) For N 0 < k„ no immune response 
can be initiated for any value of p. (C) For W 0 > k, immune response can be 
initiated for p > fS t . In this regime, the Ag-stimulation that results in 
over-suppression of immune response {/),) does not exist. (D) In this 
regime, immune response is initiated for p > p„ and is over-suppressed for 
P > p s . Note that k, and k s are dependent on the value of I. For any point in 
the plane (e.g., blue point), the values of k' s and k" are obtained by 
projecting the intersections of the line created by connecting the point to 
the origin (slope = X*) with the nonlinear curves A/ 0 = k/{X) and N 0 = k,M 
onto the W 0 -axis. By decreasing X" , the effective range of T cell activation 
(k* — k") or equivalently, the effective range of Ag-stimulation (P* s — P') 
which evoke immune response without over-suppression increases. We 
hypothesize that a healthy individual bears the potential to evoke self 
reaction and therefore its immune system is located in parametric regime 
(C) or (D); however, higher self Ag-stimulation compared to nonself 
Ag-stimulation is required for initiating immune response due to low 
renewal rate of autoreactive T cells (N 0 ). 



in the T cell repertoire (60). According to our model, a decreased 
renewal rate of a na'ive T cell clone (No) per se could prevent an 
immune response or increase the Ag-stimulation level required for 
initiation of an immune response. In addition, as shown in many 
studies, the frequency of Tregs increases with age (61, 62) which 
results in a disturbed balance between the population of na'ive T 
cells and resting Tregs (increased X). In line with these results, in 
the mathematical model an increased A. prevents the initiation of 
an immune response corresponding to the age-related immune 
hyporesponsiveness in infection and cancer. 

Based on the reasonable and physiologically realistic results that 
we could derive from the model, we dare to speculate about the 
self versus nonself concept emerging from the model. As men- 
tioned before, the na'ive T cells and resting Tregs are two major 
components of the immune reaction. The model does not dis- 
tinguish self and nonself, but rather derives different responses 
to self and nonself from quantitative differences in the nature 
of Ag-stimulation. According to the model, by adjusting differ- 
ent parameters, different requirements in terms of Ag-stimulation 



level are found for the initiation of immune responses to self versus 
nonself. If the immune system responds according to a universal 
set of Ag-stimulation thresholds, regardless of whether the stim- 
ulus arises from self or nonself-antigens, a change of systemic 
parameters can lead to immune failure or autoimmunity. Self is no 
more considered as self if it exceeds an Ag-stimulation threshold 
determined by the stringency of central and peripheral tolerances. 
Similarly, nonself is considered as self if it does not properly stimu- 
late the T cell repertoire. Autoimmunity might occur due to either 
a failure in tuning the Ag-stimulation threshold by the thymus 
that leads to unwanted self reaction in the periphery, or a chronic 
self Ag-stimulation in the periphery that leads to an oscillating 
self reaction and tissue destruction like in type 1 diabetes (63) and 
multiple sclerosis (64). Cancer or chronic infection would arise 
as the result of an imbalance in central and peripheral tolerances 
such as insufficient release of autoreactive T cells as well as high 
production or induction of Tregs that results in over-suppression 
of immune responses. 

An early elegant mathematical modeling study analyzed a series 
of models to investigate self/nonself discrimination by T cells with- 
out explicitly considering suppressive Tregs (48). As a result of their 
critical assumption that memory cells accumulate in poor stimula- 
tory conditions, the authors suggested that due to high stimulation 
by self antigens the lack of memory accumulation for T cell clones 
with high affinity to self accounts for self-tolerance. Also in our 
model, a high self Ag-stimulation (/3 > /3 S in Figure 5B) results 
in over- activation of Tregs and by this in over-suppression of self 
reaction. In both models an increased stimulation by self anti- 
gen would not lead to autoimmunity. The fact that autoreactive T 
cells do respond in the presence of Tregs when their stimulatory 
requirements are provided (15) makes it unlikely that this is the 
mechanism of self-tolerance induction. In the framework of our 
model, the view is supported that immune tolerance is induced 
by an increased stimulation threshold for self antigen and keeping 
self Ag-stimulation in a subcritical regime (jf3 < /J;). 

Undoubtedly, other mechanisms besides clonal deletion and 
Treg selection in the thymus also contribute to the fine tuning of 
the Ag-stimulation threshold required for initiation of immune 
reactions to self and nonself, such as anergy in the periphery (65) 
or activation threshold tuning in the thymus (66, 67). However, 
our simple model emphasizes the subtle balance between the gen- 
eration of Tregs and autoreactive T cells which are both needed for 
beneficial autoimmunity. The model supports the view according 
to which self and nonself do not differ on a qualitative level. It 
is rather quantitative differences of the immune status and Ag- 
stimulation level that determine which molecule is treated as self 
or nonself. 
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APPENDIX 

A.I. STEADY STATE ANALYSIS OF MODEL (1) 

The equilibrium points of model (1) can be obtained from the following: 



CTi,2i) = (Q,0) 



Equilibrium points 



(T 2 ,h) 



bf 



ad — be a 

For a positive nontrivial equilibrium point ( Ti, I2), we have to assume that: 

ad — be > 0 



(Al) 



(A2) 



The stability of equilibrium points can be determined from the sign of real part of eigenvalues of Jacobian matrix (/) . An equilibrium 
point is stable if all the eigenvalues of / evaluated at the equilibrium point have negative real parts, and it is unstable if at least one of 
the eigenvalues has a positive real part. 



Jacobian Matrix J 



al — b aT 
d-el -eT-f 



(A3) 



Characteristic Equation Q(X) = det 

i2 



X — al + b —aT 
-d+el X + eT+f 

X 2 + [eT + f + b - al] X + [-afl+ beT + bf - adT] = 0 



(A4) 



The eigenvalues (Xi^) of / for trivial equilibrium point (Ti, Ii) are obtained by solving the characteristic equation (A4): 

McTiJO : Q Wl(TiA) =^ + (f +b)X + bf = 0 



(A5) 



By checking Routh-Hurwitz stability Criterion (RHC) it can be easily confirmed that the eigenvalues have negative real parts since 
all the coefficients of polynomial Q(X) are positive, and hence, the trivial equilibrium point (T\, I\) is locally stable. For stability of 
nontrivial equilibrium point, the characteristic equation (A4) is evaluated and solved in (T2, 12): 



I 2 ( a fd \ 



X — bf = 0 



(A6) 



With the assumption (A2), the coefficient of X is positive. The sign of Q{X) coefficients change only once and hence, there exists one 
positive eigenvalue. Therefore, the nontrivial equilibrium point {T2, h) is a saddle point and unstable. 

A.2. STEADY STATE ANALYSIS OF MODEL (4) 

The equilibrium points of model (4) with definition of T cell activation (k(t)) given in equation (5) can be obtained from the following 



N - 



Np 
? + / 

dT 



bT-k 



eT+f aT 
(ad - be) T 2 + (ek - bf) T + fk = 0 
By keeping the assumption (A2), if the coefficient of T in (A9) is negative, the equilibrium points ( T), if exist, will be positive 

bf 

ek-bf<0^k<-!- 



(A7) 

(A8) 
(A9) 

(A10) 
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Otherwise, the equilibrium points will be negative. According to equation (A8), for T > 0, 1 > 0. Now, let's check the condition for 
existence of equilibrium points: 



A = 0 (ek + bf) 2 - 4afdk = (ek + bf + l/afdk} (ek + bf - l/afdk) = 0 



>o 



(* + v - 2v/*E) - o - H - ^ (. + f^)\ L. - f (> - Jl 



be 
ad 



(All) 



(A12) 



The model does not have any equilibrium points for <k <k+, and two equilibrium points otherwise. It can be verified that by 
keeping assumption (A2), we always have: 



bf 

0 <*:_< — < k+ 
e 



(A13) 



Therefore, for 0 < k < k-, condition (A10) is satisfied and the model has two positive equilibrium points and for k > k+, condition 
(A10) is not satisfied and model has two negative equilibrium points. Let's assume that the model has two positive equilibrium points 
( A > 0 and 0 < k < fc_). In the following, the linear stability of equilibrium points is analyzed: 



Jacobian Matrix J = 



P 0 



0 



p 

0 



al — b aT 
d-el -eT-f 



(A14) 



Characteristic Equation Q (X) = det 



'x + g + p 0 



0 



-P X-al+b -aT 
0 -d+el X + eT+f_ 



[X + (g + p)] [(X - a I+b)(X + eT+f) + aT (el - d)] = 0 



(A15) 



0* 



The model has one negative eigenvalue X = — (g + p ) for all equilibrium points. For the other two remaining eigenvalues, polynomial 
Q* has to be checked for existence of positive eigenvalue. 



Q* = X 2 + [eT +f + b - al] X + [-afl+ beT + bf - adT] = 0 



(A16) 



From equation (A8) it can be easily verified that b — al > 0 and hence, coefficient U is positive. Therefore, the stability depends on 
the sign of coefficient V. 



1 dT 

V = -afl+ beT +bf- adT — V 



afd 
eT + f 



T -(ad- be) T + bf 



V\ = TV ■■ 



V 2 



afd 
eT+f 



Vi 



T 2 - [(ad - be) T 2 + (ek - bf) T + fk] +k (eT + f) 



According to (A9)^=0 

dT V I= 7r+f 

+ 1 U V 2 : 



k(eT + f) kd\eT+f 
According to equations (A8) and (A9), the equilibrium values of I are: 



h = 



1 ( ( ek + bf) + ^A 

of 

1 ( (ek + bf) - 

of 



.r 2 = 



,N 2 = 



,Ni 



af , 



N 0 



b-al 2 g + p 



N 0 



b-ah g + p 



(A17) 
(A18) 

(A19) 



(A20) 
(A21) 
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where h> h and: 

A = (ek + bff - 4fkad = (ek - bff - 4fk (ad - be) (A22) 

Next, the sign of Vj, which is similar to Vi and V, has to be checked in the equilibrium points. For the larger equilibrium point 

(T2J2): 



V 2 L = ~~ ■ \e 2 k 2 + 2ekbf + ek^A + b 2 f 2 + bf^A - 4fkad] 

2fkad L J 

-sLh+^^H IA231 



V2, Vi, and V are negative for (T2, 12). Therefore, the model in this equilibrium point has a positive eigenvalue and it is locally unstable. 
For (Ti, 7i), 



V2 h 



> 0 (A24) 



V2, V\, and V are positive for (T\, h). Therefore, all the eigenvalues of the model in this equilibrium point are negative and it is 
locally stable. 

Next, let's assume that k> k + which means that the equilibrium points exist and the steady state values of T are negative, whereas 
the equilibrium values of I is positive. According to equation (A8), coefficient U in equation (A16) is negative since: 

T= — — < 0^ b- al < 0,1= dT - > 0 eT+f < 0 (A25) 

b — al eT+f 

Therefore, the model at least has one positive eigenvalue in the equilibrium points. Let's check the sign of coefficient V in the 
equilibrium points: 

V = -afl+ beT + bf - adT h ~ aI - f (al - b) + (ad - be) (A26) 

al — b 

al—b——— h 2 

V 2 = (al - b) y£=M -f(al - b) 2 + (ad - be) k = 7 - f— + (ad - be) k 

T 

T 2 (ad-be\ , 

The sign of V3 in equation (A27) which is the same as the sign of V in equation (A16) can be determined by substituting T with its 
equilibrium values from equation (A9): 

V 3 I Tl = Had [ he)fk llA + ljAjbf-ek)! (A28) 

sign? 

According to equation (A13) and reminding that k > k + , it can be verified that V3 or equivalently V is positive. Therefore, the 
equilibrium point (T2, h) has one positive eigenvalue and is unstable. For (T\, l\), 

V ^= 4(a d - b e)fk [2A + 2 ^( ek - b ^ (A29) 



V2, and equivalently, V are positive and therefore, the equilibrium point (T\, l\) has one positive eigenvalue and it is unstable. 
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A.3. STEADY STATE ANALYSIS OF MODEL (10) 

The equilibrium points of model (10) for /J = 0 can be obtained from the following 



,0, *) 
g ) 



(2i, h, NO = | 0 
N 0 

N 2 ,3 : N 2) 3 = — 



T 2)3 : ceT 2 + (cf + be - ad)T + bf = 0, 
cT 2;3 + b 



T 3 = 



T 2 = 



Ice 



Ice 



(cf + be- ad) + J~K 
(cf + be - ad) - J~A 



A = (cf + be - ad) 2 - Acfbe 



h,3 '■ h,3 = 

a 

According to equation (A32), the nontrivial equilibrium points, if available, are positive only if: 

cf + be — ad < 0 

The condition of the existence of positive equilibrium points can be obtained from equation (A32): 

A = (cf + be - ad) 2 - 4cebf = f 2 c 2 - If (ad + be) c + (ad - be) 2 > 0 
According to equation (A35), the nontrivial equilibrium points exist only if: 

{c < c_} U {c > c + } 

where 



(A30) 
(A31) 

(A32) 

(A33) 

(A34) 
(A35) 
(A36) 



{^fad — Vbej 

7 



ad - be + J**) ~ Si*) 



f 



f 



{^fad + Vbej 

7 



(A37) 



Therefore, according to conditions (A34) and (A36) and inequality (A37), the two positive equilibrium points exist only if: 

0 < c < c- (A38) 
Next, we assume that the condition (A38) is satisfied, and we analyze the stability of the equilibrium points: 



Jacobian Matrix J : 



-g-P 0 0 

P aI-b-2cT aT 
0 d-el -eT-f 



(A39) 



Characteristic Equation: 
Q(X) = det | 



x+g+p 0 0 

-P X-aI+b + 2cT -aT 
0 -d + el X + eT+f 



= & + g + P] 



X 2 + [eT + f+b - al + 2cT] X + [-afl+ 2ceT 2 + beT + bf + 2cfT- adT] 



(A40) 



All the equilibrium points have a negative eigenvalue X = —g — p. For the sign of other eigenvalues, the sign of coefficients of the 
characteristic equation has to be checked in equilibrium points. These coefficients are positive for the trivial equilibrium point and 
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therefore, all eigenvalues have negative real value. Hence, the trivial equilibrium point is stable. For stability analysis of nontrivial 
equilibrium points, the sign of U and V in equation (A40) has to be analyzed: 



U= eT + f + 2cT + (b - aI) (A33y - b a! ~ = eT+f + cT > 0 



(A41) 



Coefficient U is positive for the positive equilibrium points. Therefore, their stability depends on the sign of V in equation (A40): 

„(A33): b-aI=-cT 



V =f(b- al) + 2ceT 2 + beT + 2cfT - adT- 

= ceT 2 + [ceT 2 + (cf + be - ad) T + bf] -bf = ceT 2 - bf 



2ceT z + (be - ad)T + cfT 



According to (A32): =0 



v \r } 



1 

2ce 



A — (cf + be — ad) 



> 0 



(A42) 
(A43) 



(A44) 



According to (A34): <0 

The sign of V for the larger equilibrium point is positive and hence, this equilibrium point is stable. 



V\t 2 = 



a/A 
2ce 



\/A + (cf +be- ad) 



According to (A34): <0 



< 0 



(A45) 



According to (A32):<0 

The sign of V for the smaller equilibrium point is negative and therefore, this equilibrium point is unstable. The equilibrium points 
of the model with /S > 0 and by considering fSN(t) = k(t) are obtained by solving: 



I = 



ceT 3 - (cf + be- ad) T 2 + (ek - bf) T + kf = 0 
dT 



eT+f 



(A46a) 
(A46b) 



By keeping the assumptions (A2) and (A38), the equation (A46a) has either one positive real equilibrium point or three positive 
real equilibrium points depending on the value of k. The stability of equilibrium points that are obtained from (A46a) and (A46b) by 
varying the value of k is hard to be checked analytically; instead, it is analyzed numerically by solving equation (A40). 

A.4. STEADY STATE ANALYSIS OF MODEL (13) 

The equilibrium points of model (13) for fiN{t) = k(t) = 0 can be obtained from the following 



(Ti, h, Ru Ni) = 0, 0, 0, 



N 0 



T 2;3 : ceT 2 2 3 + (cf - ad + be)T 2>3 + bf = 0, 7 2;3 



dT : 



2,3 



, J? 2> 3 = 0, N" 2;3 



N 0 



b(eR 5 +f) b (b 

ead — be ea \e 

T 5 =0, I 5 = — ,Rs = — ,N 5 = — 
ea e g 



eT 2 , 3 + f 

e(cf + be — ad) + be(e — 1) — e(ead — be) 
cbe + y(ead — be) 



No 
g 



The sign of the equilibrium point Ti changes by changing the Treg-associated parameters (e,y). T4 is positive if 

f 

ead - be > 0,R > -- 



or 



ead — be < 0, R < — 



f 



(A47) 
(A48) 
(A49) 
(A50) 

(A51) 
(A52) 
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The physiologically relevant regime of the model occurs by satisfying the condition ( A5 1 ) which means both T4 and R4 are positive. 
The equilibrium points of model (13) for j3 N(t) = k(t) / 0 can be obtained from the following 

N=^ (A53, 
-Xk 

R = —— (A54) 
eal — b 

-l(-eXk + feaI - bf) 

(A55) 



where 



(eal - b)(-d + el) 

P 5 I 5 + P 4 I 4 + P 3 I 3 + P 2 I 2 + P\l + Po = 0 (A56) 



P 5 = a 3 fe 2 e (A57) 
P 4 = -bfe 2 a 2 e + cf 2 e 2 a 2 - ke 2 a 2 e 2 - a 3 fe 2 d - 2 a 2 fe be - a 2 e 2 X ke (A58) 

7„2, 7 _r2, .1. , „_2, i.i. , t Jr, 1.1 , .. , i.r , .i.2r. 



(A59) 

,2 f 2 



P 3 = be 2 k kea-2 cf 2 e ab + ae 2 k kb + 2 a 2 fe bd + y X kfe ae + ab 2 fe + bfe 2 a 2 d 

+ 2 ke 2 a 2 de - 2 ceX kfe a + 2 b 2 fe ae + a 2 eX ke d + 2 ke abe 2 
P 2 = 2 cek kbf -4ke abde — yk kbfe - b 2 e 2 X k - kb 2 e 2 - beX ke ad - b 3 fe - aeX kbd + cb 2 f 

+ ce 2 X 2 k 2 - ab 2 fd - ke 2 a 2 d 2 - y X 2 k 2 e 2 - 2 b 2 fe ad - y X kfe ad (A60) 
P\ = y X 2 k 2 ed + b 2 eX kd + 2 kb 2 de + 2ke abd 2 + b 3 fd + y X kbfd (A61) 
P 0 = -kb 2 d 2 (A62) 

A.5. NONLINEAR PROLIFERATION RATE OF CONVENTIONAL AND REGULATORY T CELLS 

In the models (1), (4), (10) and (13) it is assumed that proliferation rate of Tconvs andTregs is a linear function of IL-2. This simplifying 
assumption is made in order to allow parametric stability analysis of the model in a closed form and to find explicitly the dependency 
between parametric variations and topological changes of the model. Here, we show that the simplifying assumption does not affect 
the three regimes of qualitative immune responses that could be derived from the model. The linear IL-2-dependent proliferation rate 
is replaced with a nonlinear function of IL-2, named <1> (7) in models (1), (4), and (13): 



where <t> (7) is considered as a Hill-function of IL-2 



d T 

— = 0(7)7- bT + k 

dl (A63) 

— = dT-elT-fl 
dt 



r 

4>(7) = a- (A64) 

h" + 1" 



The models (1) and (A63) are compared by steady state analysis. The equilibrium points of the modified model (A63) (with 
k= g=0) are 

(T u h) = (0, 0) (A65) 
fh , / b 



^' h)= {t^' h \^h) ) (A66) 

The nontrivial equilibrium point (T2, h) is positive and biologically meaningful only if 

(a- b)d n - be n h n > 0 (A67) 
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The local stability of the equilibrium points can be determined by obtaining the eigenvalues from the characteristic equation: 
Characteristic Equation Q(X) = det 



X-aJ^ + b -aT nh " in 



X 2 + 



-d+ el 
eT+f - a 



jn 



h" + I" 



(h"+I"f 
X + eT+f 



+ b 



(A68) 



+ 



\- a T^TTn + b ) ( eT + f) + aT V, 772 (- d + eI A=° (A69) 

V h" + 1" ) v ' (h n + I») 2 J 



By checking Routh-Hurwitz stability Criterion (RHC) it can be easily confirmed that the eigenvalues have negative real parts for 
trivial equilibrium point since all the coefficients of polynomial Q(X) are positive, and hence, the trivial equilibrium point (T\, I\) is 
locally stable. For checking the stability of the nontrivial equilibrium point, the characteristic equation (A69) is evaluated in (T 2 , I 2 ): 



Xi,2\(T 2 ,i 2 ) ■ QW\(T 2 ,I 2 ) = X 2 + [eT 2 +f] X + 



nbh n T 7 



h(h" + I 2 n ) 



(d - eh) 



0 



(A70) 



By assuming the condition (A67), the coefficients U and V are positive and negative respectively. Therefore, the sign of the coef- 
ficients of Q(X) (U and V) changes only once and hence, there exists an eigenvalue with positive real part. Therefore, the nontrivial 
equilibrium point (T 2 , I 2 ) is a saddle node and unstable. 

Similar to the model (1), the stable manifold of saddle node in the model (A63) defines a threshold for the initial conditions that 
allow for unlimited proliferation of activated T cells. By comparing the conditions (A67) and (3), the dependencies of these conditions 
to the model parameters, specifically the proliferation rate (a) and IL-2 secretion rate (d), are positively correlated. In other words, 
in both models, only T cell clones with sufficiently high proliferation rate (a) and/or high IL-2 secretion rate (d) are able to undergo 
major T cell proliferation. 

For k / 0, the equilibrium points of the model (A63) are obtained from the following equations: 



fl 



d- el 



I: [(a- b)f - ke] + [kd] I" - [h n (bf + ke)] I + kdh" = 0 



(A71) 
(A72) 



The stability of equilibrium points is analyzed by evaluating the characteristic equation and is shown in Figure Al for parameter 
values given in Table 1 and Hill-function parameters n = 2 and h = 0.5. By comparing the bifurcation diagram in Figures Al and 3B, 
the qualitative similarity between model (4) and (A63) is evident. This qualitative similarity also holds true between model (13) and 
the following model: 



dT 

d7 

dR 

df 

d/ 

dt 



= <D(J)T - bT - cT 2 - yRT + k 
= e<t>(I)R- bR + fiN 
= dT - eI(T + R) - fl 



(A73) 



where 4> (I) is identical to (A64). 

The equilibrium points of model (A73) are calculated and their stability is analyzed by deriving the characteristic equation of the 
model and obtaining the eigenvalues. By keeping assumption (12), the bifurcation diagrams of model (A73) for two different values of 
X are obtained by treating k = > 0 as the bifurcation parameter (depicted in Figure A2). Depending on the value of k, the model 
has either 8 or 6 equilibrium points (4 or 2 equilibrium points with T > 0, identical to model (13)) with parameter values given in 
Table 1 and Hill-function parameters n = 4 and h=l. The additional equilibrium points resulted from considering the Hill- function 
nonlinearity are all in the negative space of the model variables. As it can be seen from Figure A2, similar to the model (13), the three 
qualitatively different responses still could be derived from the modified model. It is clear that the value of fc;, k s , and X t h are different 
from their corresponding values in the model (13). 

In summary, imposing the nonlinear IL-2 dependent proliferation rate of cells results in a more restricted condition for initiation 
of an immune response in comparison to the linear IL-2 dependent proliferation rate, namely the requirement of higher T cell avidity 
(higher a and d), higher Ag- stimulation (increased ft ,), and lower Treg/Tconv ratio (lower A. t j,); but three qualitatively different immune 
reactions depending on the critical levels of Ag-stimulation could still be derived, very similar to model (13). 
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10 2 




Steady state value of T cell activation k 



flN„ 
g + P 



FIGURE Al | Bifurcation diagram of model (A63) with Hill-function parameters n = 2 and h = 0.5 by treating k as bifurcation parameter Stable and 
unstable equilibrium points are shown by black and red lines, respectively. For k > the immune response enters the regime of unlimited proliferation. The 
unstable negative equilibrium point is omitted in this figure. 
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FIGURE A2 | Bifurcation diagram of model (A73) with Hill-function 
parameters n = 4 and h = 1 using k as the bifurcation parameter with 
(A) I = 0.0016 and (B) A. = 0.0008. Stable and unstable equilibrium points are 
drawn by black and red solid lines, respectively. The stable limit cycles are not 
shown for all value of k f <k< k s except for k= 12. Depending on the values 
of X and k, an immune response is not initiated (gray), is initiated (red) or 



over-suppressed (magenta). With parameter values given r Table 1 and 
Hill-function parameters n = 4 and h= 1, the threshold becomes A th = 0.00111. 
The time-courses of the activated T cell population 7~(f) were deduced from a 
numerical solution of model (A73) with zero initial conditions and persistent p. 
The negative equilibrium points which are all unstable are not shown in the 
plots. 
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